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Abstract 

A method for making realistic estimates of the density of levels in even-even 
nuclei is presented making use of the Monte Carlo shell model (MCSM). The 
procedure follows three basic steps: (1) computation of the thermal energy 
with the MCSM, (2) evaluation of the partition function by integrating the 
thermal energy, and (3) evaluating the level density by performing the in- 
verse Laplace transform of the partition function using Maximum Entropy 
reconstruction techniques. It is found that results obtained with schematic 
interactions, which do not have a sign problem in the MCSM, compare well 
with realistic shell-model interactions provided an important isospin depen- 
dence is accounted for. 
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The density of levels in nuclei plays an important role in understanding compound nuclear 
reactions. Two particularly important examples are the decay of the giant-dipole resonance 
(GDR) in hot nuclei and the radiative capture of light nuclei, i.e., protons, neutrons, 
and alphas, in nucleosynthesis [[J . In the first case, properties of the GDR, in particular the 
damping width, have been studied in several nuclei for excitation energies ranging from 50 to 
200 MeV, and it has been shown that the analysis of experimental data is very sensitive to the 
the dependence of the level density on excitation energy |§. In contrast to the GDR studies, 
the particle capture probability, which determines the rate at which nucleosynthesis reactions 
occur, is sensitive to the level density near the particle-decay threshold: i.e., ~ 5 — 15 MeV. 

In most applications where the level density is required, the Fermi-gas model estimate 0] 
is employed 

where E is the excitation energy, and a is the level-density parameter, which is determined 
by the number of single-particle states at the Fermi energy. The principal shortcoming of 
the Fermi gas estimate is that interactions between nucleons are ignored. Effects due to shell 
corrections and pairing correlations are approximated in Eq.([TJ) by replacing the excitation 
energy E by backshifted quantity £-A|. Emprically, both a and A exhibit a dependence 
on E and the number of nucleons, A, that cannot simply be estimated within the context 
of the Fermi-gas model; a typical value for a at low excitation energies is a ~ A/8. 

An alternative model that explicitly includes both one- and two-body correlations is the 
shell model. State-of-the-art shell-model Hamiltonians, such as the universal SG?-shell (USD) 
Hamiltonian of Wildenthal have been very successful at describing both excitation en- 
ergies and transition amplitudes for states in a wide range of nuclei (18 < A < 48) up to 
excitation energies of the order 5-10 MeV. Although the shell model might appear to be the 
obvious method for estimating the level density at low excitation energies, direct diagonal- 
ization of the Hamiltonian faces severe computational limitations due to the fact that the 
number of basis states scales as the exponential of the number of valence particles. Indeed, 
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the large number of basis states was the motivation behind the development of the spectral 
distribution methods of French and others JF], which rely on the statistical properties of the 
elements of the Hamiltonian matrix. A central-limit theorem can be applied to describe the 
action of the Hamiltonian in large spaces ||; generally reducing the problem to the calcu- 
lation of the first and second moments of the Hamiltonian. In many practical applications, 
however, this limiting situation may not be sufficiently realized. For example, because of 
features of the Hamiltonian, which may be thought of as shell corrections, significant depar- 
tures from "normality" may be observed at low excitation energies. As a consequence, it is 
necessary to compute higher-order moments of the Hamiltonain or partition the shell-model 
space into smaller subspaces. Unfortunately, not only are these higher-order moments more 
difficult to evaluate, but the level density reconstructed with orthogonal Hermite polyno- 
mials may fail to be positive definite ||. An alternative method was proposed by Pluhaf 
and Weidenmiiller |L0| in which the partial level densities in the subspaces were assumed 



to have a form predicted by the Gaussian orthogonal ensemble (GOE), i.e., semicircular, 
as determined from the first and second moments of the Hamiltonian within the projected 
subspaces. The total level density, which in their method is guaranteed to be positive defi- 
nite, is then obtained by combining the various subpartitions with the coupling between the 
subspaces being determined statisitcally from the mean off-diagonal matrix elements. This 
procedure also faces several limitations because of the reliance on the GOE limit for the 
subspaces, as well as being restricted to the first and second moments of the Hamiltonian. 
Consequently, in applications to more general shell-model problems, this procedure tends to 
lead to level densities that are somewhat broader than the exact results [|ll[ . 



In this work, a method for making realistic estimates of shell-model level densities using 
the Monte Carlo shell model fl2| (MCSM) is presented. The power of the MCSM is that 
it is capable of providing exact results for a range of observables in model spaces where the 
dimensions are prohibitive for direct diagonalization. In addition, the MCSM is quite well 
suited to compute thermal properties, such as the energy, from which the partition function 
may be obtained, which then yields the level density through an inverse Laplace transform. 



The applicability of the MCSM to the most general of shell-model Hamiltonians, however, is 
limited because of the sign problem associated with the Monte Carlo weight function. One 



is then faced with either using an extrapolation method |Tj|, which tends to yield larger 
statistical errors, or a schematic interaction that is free of the sign problem. Here, it will be 
shown that schematic Hamiltonians, such as a surface-delta interaction |TJ| , possess most of 



the global, or collective, features exhibited by "realistic" Hamiltonians with one important 
exception: in the spectra of even-even nuclei, the higher isospin states tend to be too low 
in energy, thereby compressing the total level density. This improper isospin dependence 
can be corrected by adding the term aT 2 to the schematic interaction, thus shifting the 
excitation energy of the higher isospin states. Unfortunately, this T 2 term also has a bad 
sign, and cannot be computed directly. To address this problem, a simple and accurate 
approximation for correcting the thermal energy for the T 2 dependence in even-even nuclei 
is presented, and it will be shown by direct comparison that an MCSM calculation using a 
schematic interaction yields a reasonable estimate of the level density obtained in "realistic" 
shell-model calculations. 

The procedure consists of four steps: (1) using a semi-realistic schematic interaction, 
compute the thermal expectation value of the Hamiltonian, 

Tr[He-W] E i E l e~^ fdEe^ E Ep(E) 
[P) Tr[e-^] E»e-^ fdEe-? E p(E) ' { ' 

with the MCSM; (2) correct E(/3) for the missing T 2 dependence using Eq. (^) below; (3) 
compute the partition function, Z(/3), via 

In Z(B) = - I dp'E(p') + In Z(0), (3) 
Jo 

where Z(0) is actually the total number of states; and (4) calculate p(E) using maximum- 
entropy reconstruction techniques to perform the inverse Laplace transform of Z((3). 

The primary goal of this work is to establish the feasibility of the method for making 
estimates of the level density in realistic situations. To accomplish this, comparisons with 
exact results are necessary, and the focus of this work is the nucleus 24 Mg, which has four 



valence protons and neutrons occupying the 0d 5 / 2 , lsi/2> and 0d 3 / 2 valence orbitals. Because 
of the overall success of the USD interaction, which was explicitly developed for this model 
space, this interaction will be used as a benchmark for success. The schematic interaction 
is composed of the three single-particle energies and a two-body potential given by the 
surface-delta interaction (SDI) [TJ[], which has the form 

U(n, r 2 ) = -47ry 5(r 1 - r a )«J(n - Ro), (4) 

where r« is the position vector for the i th particle and Rq is the nuclear radius. The prin- 
cipal feature of the SDI is that it is basically comprised of multipole-multipole terms, 
with the dominant multipoles being monopole and quadrupole. Taking R = 3.145 fm, 
Vq = 54.76 MeV/fm 2 , and evaluating the two-body matrix elements using harmonic oscilla- 
tor wave functions with hu = 13.531 MeV, the T = 1 USD matrix elements are reproduced 
by better than 500 keV. The single-particle energies are adjusted to reproduce the low-lying 
experimental spectra for 19 and the USD shell-model spectra for 25 and 27 0. Using fixed 
single-particle energies across the shell and the USD mass scaling of (18/A) 3 for the two- 
body part, the SDI Hamiltonian reproduces the USD spectra and level densities for oxygen 
isotopes reasonably well. On the other hand, an important component is missing, as it is not 
possible to reproduce the binding energies for the oxygen isotopes. This feature is generic 
to schematic interactions of the form multipole plus pairing, and may be "fixed" by adding 
a term dependent on the square of the isospin, aT 2 , as is illustrated in Fig. [I], where the 
binding energies for the sd-shell oxygen isotopes are compared for the USD interaction and 
the schematic Hamiltonian, with parameters = —4.820 MeV, e s i/ 2 = —2.820 MeV, 
e d3/2 = 1.530 MeV, 1/ =54.76 MeV/fm 2 , and a = 0.546 MeV. 

From the standpoint of estimating the level density, an adequate Hamiltonian would in 
fact be an extension of the SDI that includes isospin-dependent components: in particular 



the modified surface delta interaction |T5| . Because of the sign problem in the MCSM, 
however, it is generally not possible to treat isospin-dependent terms in the Hamiltonian. 
On the other hand, for a given isospin value, the SDI spectra compare well with the realistic 



USD results. Hence, a reasonable first-order correction may be obtained by simply shifting 
the higher isospin states by adding aT 2 to the Hamiltonian. The magnitude of the coefficient 
a, however, cannot be estimated a priori, as it seems to differ from nucleus to nucleus and 
depends on the number of the protons and neutrons as well as whether the shell is more 
than half full. This is illustrated in Table |IJ, where the excitation energies obtained with 
the SDI and USD interactions for the first J = + states of each isospin are tabulated for 
24 Mg ans 32 S. Also, given in the table are the total number of states of each isospin that 
contribute to the total level density. Because of this unpredictable T 2 behavior, a shell- 
model calculation with a realistic interaction, even within a truncated model space, for a 
single angular momentum for a few isospins is necessary to determine a reasonable estimate 
for a. For 24 Mg, the improper isospin dependence is corrected on average with a = 1.507. 
On the other hand, to illustrate that the SDI interaction gives a good representation of the 
low-lying collective behavior for 24 Mg, the excitation energies of the lowest few T = states 
are compared with the USD values in Table |TJ. 

Using the Monte Carlo shell model (MCSM) techniques described in Ref. [12], the ex- 



pectation values of observables such as the Hamiltonian, T 2 , etc., were evaluated for the 
SDI interaction as a function of the inverse temperature, (3, in the range < (3 < 1 MeV -1 
in increments of A(3 = 1/80 MeV -1 . In order to ensure sufficient accuracy, 2000 Monte 
Carlo samples were taken for each (3 value, and multiple time slices were used in all MCSM 
calculations with (3 > 0.0375 MeV -1 , with the maximum number of forty being used at 
(3 = 1 MeV -1 . Typical Monte Carlo uncertainties for the energy observable ranged from 
10 keV for small values of (3 to about 120 keV for /3 ~ 1 MeV -1 . 

With E((3), the partition function is obtained via Eq. @, and the level density is then 
given by the inverse Laplace transform of Z((3). Although a saddle-point approximation 



may be employed, giving p(E) = e^ E+lnZ /J— 2irdE / d(3, this method tends to be somewhat 
unstable at low excitations energies due to difficulties associated with computing the deriva- 
tive in the denominator. An alternative method is to evaluate the inverse transform using 
maximum-entropy (MaxEnt) reconstruction techniques [16|]. The starting point is to bin 
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p(E) into Nr bins of equal width AE, namely 

N R 

p(E) =J2fi [8( E ~ ( E i ~ A /2)) " 0(E - (E { + A/2))] , (5) 

i 

where fi is the number of levels contained within the i bin. With this level density, the 
reconstructed partition function is then 

N R O f 

Zr(P) = E ~jr^ Ei sinh(/3A£/2). (6) 

The goal of MaxEnt is to find the set of values {fi} that maximize the extended entropy 
functional aS — x 2 /2, where 

j °i 

quantifies how well the reconstruction reproduces the N calculated values of the Z(/3j) (note 
N r < N) and the information entropy, S, is given by 

Nr 

S = J2(fi- D i~fiHfi/Di))- (8) 

i 

In the MaxEnt method, it is necessary to specify a default model {Dj}, which may be used to 
characterize any prior information known about the problem at hand. In this case, it is well 
known that within a finite model space, the level density exhibits a Gaussian character [ITT 



which may be used to define the default model. For finite-space, shell-model calculations, the 
total number of states is known and the first and second moments of the level density may 
be obtained from E((3 = 0) and dE((3)/dP\p=o, respectively. The reconstructed {fi} also 
depend on a, which governs the relative weight between the default model and chi-square. 
Here, a was chosen so that x 2 ~ N. Finally, the uncertainty in the fi values may be obtained 
in a manner similar to least-squares fitting from the curvature matrix d 2 (aS — x 2 : /2)<9 ] fid fj. 

As was mentioned above, the SDI interaction exhibits a bad isospin dependence that can 
be corrected by adding the term aT 2 to the Hamiltonian. Unfortunately, this additional 
term has a bad sign in the MCSM and cannot be evaluated directly. The extrapolation 
method of Ref. [TBI] could be used, but at a significant computational cost and larger sta- 



tistical errors. Instead, for even-even nuclei, where the low-lying states are unaffected by 



the additional term, a perturbative approach may be more useful. Towards this end, aT 2 is 
added to the SDI Hamiltonian, and the energy is evaluated by expanding the aT 2 terms in 
the exponential in Eq. (Q) to first order in (3. Unfortunately, a further limitation is imposed 
due to computational limitations that make it impractical to evaluate the expectation value 
of n-body operators in the MCSM beyond n = 2. Given these considerations, the first-order 
correction to the energy is estimated as 

E corr = (3a^^ + a(f 2 ). (9) 

In comparison with exact results, Eq. (H) works quite well, although it tends to "over correct" 
by approximately 10%. This over correction, which may be due to the neglected (T 2 ) 2 terms 
in the expansion, can be damped by multiplying E corr by the factor e~ /3a ^ T2 K For illustrative 
purposes, both E{(3) and In Z(f3) are shown in Fig. |2] for 20 Ne using the SDI interaction with 
and without the iosospin correction factor. In the figure, the solid and dashed lines represent 
E(/3) and Z(f3) obtained by using the exact eigenvalues of a shell-model diagonalization for 
the SDI and SDI + at 2 Hamiltonians, respectively. The dotted line represents the results 
obtained by adding E corr to the SDI results. From the figure, it is seen that the corrected 
values very accurately reproduce both the energy and partition function of the SDI + aT 2 
Hamiltonian. 

Shown in Fig. ^|are the reconstructed level densities obtained for 24 Mg using the SDI (top) 
and SDI + 1.507T 2 (bottom) Hamiltonians, respectively. For the reconstruction, E(0) was 
normalized relative to the ground state value, which was evaluated to be —76.844(200) MeV 
by performing MCSM calculations for 1.75 < (5 < 2.25 MeV -1 and extrapolating for f3 — > 
oo. The histogram in both panels represents the exact level density for 24 Mg using the 
USD Hamiltonian, and was obtained by direct diagonalization and placing the 28503 (J z = 
0) eigenvalues into 1 MeV bins while including the 2J + 1 degeneracy. Because of the 
near Gaussian structure of the level density, p(E) is only plotted up to 45 MeV, which 
is slightly larger than the centroid of the USD level density (~ 42 MeV). From the top 
panel, it is clear that the level density obtained from the SDI Hamiltonian alone is much too 
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compressed (centroid at ~ 36 MeV), and over predicts the p(E) by nearly a factor of two 
for E < 20 MeV. On the other hand, the corrected level density represents a considerable 
improvement (centroid ~ 40 MeV), and gives a reasonable representation of the USD level 
density for excitation energies < 20 MeV. 

At this point, it is important to note an important limitation in using the shell model to 
estimate the nuclear level density. Because of the fact that all calculations are by necessity 
limited to a finite model space, the shell-model level density will always have a Gaussian 
character, and at some point will always underestimate the true level density because of the 
presence of states representing excitations outside of the model space. For the most part, 
the shell model is best suited to describe states at lower excitation energies. Hence, it is 
not realistic to expect that the shell model could provide an estimate of the level density for 
excitation energies around 40-50 MeV. In the case of a single major oscillator shell, such as 
the sd-shell used here, the results are comparable to experimental data for E < 10 MeV. On 
the other hand, by enlarging the model space to include more configurations, say another 
major oscillator shell such as the /p-shell, the method presented here can be used to make 
estimates for the level density up to excitation energies of the order 15-20 MeV. In this case, 
however, it will also be necessary to account for states of differing parity, as well as spurious 
excitations of the center of mass. 

To conclude, the procedure outlined here can successfully describe the level density of 
a realistic nuclear system, and plans are currently underway to compute level densities 
for 22 Mg, 26 Si, 30 S, and 34 Ar, which are needed to make estimates of (a, 7) reaction rates 
of astrophysical interest ||18|| . In addition, the method will also be applied to 162 Dy and 



172 Yb where experimental data for p(E) exists for excitation energies up to the neutron 



separation energy. In these nuclei, since the protons and neutrons occupy different major 
shells, the problems associated with the higher isospin states pointed out here are most likely 
to be mitigated. 

Discussions with C. W. Johnson and J. P. Draayer are gratefully acknowledged. This 
work was supported in part by NSF Cooperative agreement No. EPS 9550481, NSF Grant 
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TABLES 

TABLE I. Excitation energies (in MeV) for the first J = states of each isospin for the SDI 
and USD interactions for 24 Mg and 32 S. The number of states for each isospin contributing to the 
total level density is also tabulated. 



T 


24 Mg 

USD 


SDI 


USD 


32 S 


SDI 


# States 





0.000 


0.000 


0.000 




0.000 


166,320 


1 


12.872 


9.351 


7.312 




4.735 


332,640 


2 


15.425 


7.903 


12.060 




6.434 


237,600 


3 


33.923 


22.528 


33.468 




22.268 


83,160 


4 


46.040 


26.251 


45.404 




26.328 


11,880 



TABLE II. Comparison between the USD and SDI excitation energies (in MeV) for the lowest 
few T = states in 24 Mg. Shown are the first ten states as predicted by the USD interaction, and 
then the lowest J = 7 and 8 states. 



J USD SDI J USD SDI 

0.000 0.000 7.561 6.778 

2 1.509 1.973 1 7.764 7.253 

2 4.122 4.146 5 7.883 8.281 
4 4.378 4.845 6 8.263 8.460 

3 5.097 5.326 7 12.283 12.596 

4 5.934 5.314 8 12.088 10.974 
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FIGURES 

FIG. 1. Comparison of binding energies for oxygen isotopes as a function of neutron number 
obtained with the USD (solid line) and SDI + 0.546T 2 (dotted line) Hamiltonians, respectively. 

FIG. 2. Plot of the energy, E(/3), and partition function, Z{(3) for 20 Ne, as a function of (3 for 
the SDI and SDI + aT 2 Hamiltonians. The solid and dashed lines represent the SDI and SDI + aT 2 
results obtained from the exact eigenvalues, while the dotted line shows the SDI results corrected 
via Eq. (g). 

FIG. 3. Reconstructed level densities obtained for 24 Mg using the SDI (top panel) and the 
SDI + 1.507T 2 (bottom panel) Hamiltonians, respectively. The histogram in both panels represents 
the 24 Mg level density obtained with the USD interaction. 
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